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, The model Hamiltonian of a two-dimensional Bose liquid (proposed earlier by Kane, Kivelson, Lee 

' and Zhang as the Hamiltonian which has Jastrow-type wavefunctions 'I'o{ri} oc Jl^^;, l^j — rfc|^" 

as the ground-state solution), is shown to possess nonrelativistic supersymmetry. For the special 
» I ' value of the coupling constant a = 1/2 the quantum mechanics described by this Hamiltonian is 

' shown to be equivalent to the dynamics of (complex) eigenvalues of random Gaussian ensemble of 

normal complex matrices. For general a, an exact relation between the equal-time current-current 
and density-density correlation functions is obtained, and used to derive an asymptotically exact 
lO I (at low wavevectors q) spectrum of single-particle excitations beyond the superfluid ground-state 

. (realized at low a's). The ground-state at very large a is shown to be of "Quantum Hexatic" 

type, possessing long-range orientational order and quasi-long-range translational order but with 
zero shear modulus. Possible scenaria of the ground-state phase transitions as function of ct are 
■ discussed. 



PACS numbers: 05.30.Jp, 71.10.Hf, 12.60.Jv. 



I. INTRODUCTION 



It is commonly believed that the system of Bose particles at zero temperature may exist only in one of two possible 
' ground-states: superfluid (SF) or crystalline (CR). Most probably it is indeed the case if the interaction between 
bosons is short-range (compared to the mean interparticle distance n~^/^). However the situation is much less clear if 
O long-range and/or long-time interactions are involved. Two physical situations where the problem of a two-dimensional 
(2D) Bose system interacting with a 2D Maxwell gauge field were considered in ||l|: classical thermodynamics of Vortex 
Line Liquid in High- Temperature Superconductors (see for an extensive review) and low-temperature quantum 
properties of liquid of Holons (arising in the RVB-type theories [|| of doped Mott dielectric Q). The arguments were 
given [Ql in favor of existence of a new unusual ground-state which is still liquid, but is not superfluid. More recently 
the same problem was attacked by another approach also in relation with RVB-type theories. Still the theoretical 
^vq ' arguments in favor of the existence of such unusual ground-state are not quite conclusive; the results of different 
I computer simulations are in mutual disagreement [^|j7j. Whereas in two phase transitions and an intermediate 
' state were observed, in Q a direct first-order transition from SF to CR state was found; note however, that the model 
Hamiltonians simulated in these two numerical studies were different in that the retarded interactions were presented 
5!!;; onlyin§l. 

. There are two different questions related with this problem: i) can one invent such a ground-state in a self-consistent 
way? and ii) is it realized in some model with a physically realistic long-range interaction? In the present paper we 
I try to approach the answer to the first of the above questions, using a special model of a 2D Bose liquid, introduced 
originally by Kane, Kivelson, Lee and Zhang (KKLZ) The KKLZ model Hamiltonian is constructed in such a 
C way that its exact ground-state wavefunction is of the Jastrow type, i.e. it is a product of 1-particle and 2-particle 
^ ' dependent factors. With the special choice of the Jastrow function the main (in the long-wavelength limit) term of 
• • , the KKLZ Hamiltonian is given by a pair-wise logarithmic interaction, i.e. it resembles a 2D static Coulomb potential 
of the real problem discussed in [|l|. On the other hand, special symmetry (in fact, a kind of a supersymmetry) of the 
KKLZ Hamiltonian makes it possible to obtain a number of exact relations between the equal-time T = correlation 
$H functions, which help to analyze the general structure of the theory. 

The paper is organized as follows: in Section II we derive the KKLZ Hamiltonian from the condition that it is 
given by the zero-Fcrmion sector of a supersymmetric nonrelativistic parity-conserving theory, which resembles (in 
the long-wavelength limit) the system of 2D bosons with logarithmic interaction. The coupling constant a of our 
theory is proportional to the charge of 2D bosons, thus the ground-state possesses superfluidity at low a . In Sec. Ill 
we analyze a special case (a — 1/2) of the obtained model and show that it is equivalent to the quantum mechanics 
of a Gaussian Matrix Model with normal complex matrices (in the same sense as quantum mechanics of Gaussian 
Hermitean Matrix Model is equivalent to ID free fermions |^); in Sec. IV we are back to the general case and derive an 
exact real-space relation between the equal-time current-current and density-density correlation functions (extending 
the similar results found in for the Calogero-Sutherland model to a 2D case). These relations are then used in 
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Sec. V in order to derive a relation between momentum-space transverse current-current correlation function and the 
static structure factor, and in Sec. VI to find an exact form of the low-momentum excitation spectrum u>{q) (which 
appears to be given by the one- loop correction to the MFA). Sec. VII is devoted to the discussion of the "crystalline" 
ground-state of the KKLZ Hamiltonian; we show that such a state has in fact zero shear modulus and thus should be 
considered as a kind of a "quantum hexatic" (QHX), in some analogy with the classical hexatic known to exist as an 
intermediate phase within the Halperin-Nelson picture of thermal 2D melting jl^. In Section VIII we speculate on 
possible properties of the additional ground-state which may be sandwiched between the SF and QHX states. Our 
conclusions are presented in Sec. IX. Finally, Appendix is devoted to the description a model of classical stochastic 
dynamics which is equivalent to our quantum-mechanical problem; we show how to use this equivalence for efficient 
numerical simulations of some quantum-mechanical expectation values. 



II. THE KKLZ HAMILTONIAN: SUSY STRUCTURE AND GENERAL FRAMEWORK 

We will see now that the KKLZ Hamiltonian can be considered as an example of many-body generalization of the 
supersymmetric quantum mechanics introduced in |l^ , [l3|| . We define a pair of the conjugated "supercharge" operators 
Q and Q+i 

Q = Y^ Iho.O'tc'^ = H 9j>ai,a; (1) 

where Fermionic operators aj^a and ^ corresponds to each component of the j-th particle's coordinate vector rj^a 
(with a = 1, 2), and 



. 9 „ /v-^ 91n Ir," — rfcl \ 



Note that the introduction of two Fermionic operators aji, aj2 is necessary in order to construct a P-invariant theory; if 
the condition of P-invariance were relaxed, the use of just one (per each particle) Fermionic operator Cj corresponding 
to the holomorfic coordinate zj = xj + iyj would be possible. Now we define the Hamiltonian 

^ = ^W,Q+} (3) 

and note that this Hamiltonian, when acting on the state of "Fermi-vacuum" (defined by aj alOp) — 0), H = T-C\Of), 
is equivalent to the Kane-Kivelson-Lee-Zhang (KKLZ) Hamiltonian |^ (see also [|l|) 

H~E° = Fkklz = ^ + V2 In |r, - r, | + 4a^ ^ V, In |r, - r, |V, In |r, ~ rfe| + {2TTanfA, 

(4) 

where — —2TranN[l + a{N — l)]/m. The corresponding ground-state wavefunction is given by 

Mri, tat) = const • J] - rfel^^e^ 'E.-^ (5) 
j>k 

Alternatively, the Hamiltonian can be written directly in terms of the qj^a operators: 

j 

Each of the Qj^a operators annihilates the ground state (^, so the ground-state energy corresponding to H is zero. The 
Hamiltonian (^) contains 2-body and 3-body interactions; in terms of the local particle's density p(r) the corresponding 
potential energy is given by 
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U=IJ d^rid^r2p{ri)p{r2)V2iri-r2) + J d^rid^r2d^r3p{ri)p{r2)pir3)V3{ri - r2,ri ~ r^), (7) 



where 



V^2(r)- ^(r) + l^3(r,r') = 

Writing p(ri) = rt + 8p(v\) where n is the niean density, and integrating over ri in the three-body term of Eq. we 
an find effective two-body potential 

2 

Vf{r)^-^\r^v^V2{r\ (8) 

ZTT 



where we introduced the Coulomb coupling constant g — ^■na\J njm in the same way as in 0. Now the potential 
energy can be rewritten in terms of the density variations 5p — p ^ n, with the replacements p(r) — > 5p{v) and 
V2{y) Vf^(Y) made in Eq. (0). On large space scales the main term in the potential energy is the one containing 
Inr, which leads to the existence of a plasmon gap in the longitudinal branch of the excitation spectrum: uj{q = 0) = 

Below we will be interested in the presence of superfluidity in the ground state defined in Eq. (^. We define 
superfluid density via the variation of the ground-state energy Eq (do not be confused by the different quantities 
Eq and in Eq. (^) with respect to the external (source) gauge field a minimally coupled to the Hamiltonian (^) 
via the replacement Vj Vj + ia(rj) in the kinetic energy term only: 



TO aaa(q)5a^(-q) 



■ {Sc.f3 - q»qf3) ^ riq^O,u; = 0), (9) 

q~*0 



where the polarization function V{q,u} — 0) is related to the transverse current-current correlation function Ii±{q,uj) 
by the zero-T FDT relation: 

nq,. = 0)^'-r^^^d^. (10) 

71" Jo ^ 

Below in Sec. IV we will find the relation between the instantaneous transverse current-current correlation function 
n_L((z) = ^ /g°° nj_((7, Lj)duj and the density-density correlation function. Note that the replacement V — > V -I- ia made 
in all terms in Eq. (|^) would retain SUSY structure of the theory and therefore would not change the ground-state 
energy Eq = 0. 

The ground-state wavefunction is given by Eq. (^) for any value of the coupling constant a, however this multiparticle 
wavefunction describes at least two ground-states which are very different physically. Namely, at a < 1 the ground- 
state is a superfluid |^ and its superfluid density is equal to the total particle's density n. However, superfluidity 
is not guaranteed by the form of the ground-state wavefunction only: at very large values of a (larger than acr ~ 35) 
the ground state is expected to possess translational symmetry breaking realized in the formation of a triangular 
lattice (in fact only quasi-long Berezinsky-type |Q translational order is formed). Such a state does not possess 
superfluidity, Us = 0. Therefore at least one quantum phase transition, between the states with nonzero and zero 
Us, should happen as the parameter a increases. There is also a possible room for the second phase transition. To 
see it, note that the Galilei invariance of the Hamiltonian ^ implies (at T = 0) the absence of a normal density: 
Un = n — Us = 0. It means that nonzero n„ can appear due to spontaneous breakdown of the Galilei invariance, 
which may happen simultaneously with the lattice formation (so that Ug jumps from n to zero), or at some other 
critical value ac < ctcr- Moreover, there are no special reason for the superfluid density to go to zero at a = acr, as it 
may happen at a G [ac, ctcr]- Finally, there is one more symmetry, the orientational one, which can be spontaneously 
broken by the ground-state at a = ao G [aciCkcr]- Below we start our consideration from the small-a range where 
no symmetry is broken, then proceed (in Sec. VII) to the region of largest a > acr, and at the end we present some 
speculations on the system's behavior in the intermediate range of a. 

III. THE SPECIAL CASE a = 1/2: MAPPING TO QUANTUM MECHANICS OF NORMAL MATRICES 

We will show now that the dynamics of a Gaussian ensemble of iV x TV complex matrices {N — > oo) constrained 
by the condition [M, M+] = (which is the definition of a normal matrix) is intimately related to the many-body 
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quantum mechanics defined by the Hamiltonian (Q) with a = 1/2. This relation is of the same origin as the well-known 
now relation between dynamics of the Dyson matrix ensembles and the Calogero-Sutherland model of many-body 
quantum mechanics in 1 dimension, where the coordinates of ID particles are just eigenvalues of the corresponding 
matrices. Eigenvalues of a normal matrix are, in general, complex numbers, which can be considered as coordinates 
of particles living in a 2D space. Thus our goal here is to show how it is possible to perform functional integral over 
most part of variables characterizing dynamics of a normal matrix and reduce the problem to the quantum mechanics 
for the set of its complex eigenvalues Zi^Zi. Our derivation follows closely the one presented in (see also for 
three Dyson ensembles. Consider the partition function 



Z = I PM(r)PAf+(r)exp[--S'(M,M+)][|(5([M(T),Af+(r)]) 
with the flat measure of integration and the Gaussian action 



5(M,M+) 



1 dM+ dM 

- Tr 

4 \ dr dr 



+ Tr (M+M) 



dr. 



(11) 



(12) 



Note now that any normal matrix can be diagonalized by an appropriate unitary transformation: M = U^ZU , where 
Z is diagonal with the elements z^. Therefore we can change variables in the functional integral (11): 



VM{t)VM+{t) = J{z,}Y[dz,dziVU, 



(13) 



where the Jacobian of the transformation is given ||lf 
determinant: J{zi} = const |A{2:i}p = const Yl^-^j \zi 
above functional integral in the transformed variables: 



by the square of the absolute value of the Vandermonde 
Zjp. Consider now the Hamiltonian corresponding to the 



d d d d 

dzi dzi dzi dzi 



El 



H. 



(14) 



Here the first sum over i constitutes a Laplacian operator in the curved space with the square root of the determinant 
of the metric tensor y/g — |A{zi}p, the second one comes from the "potential energy" Tr (M+M), and the last term 
corresponds to rotations of angular (i.e. related to unitary transformations) variables, and thus vanishes when acting 
on any singlet (J7-invariant) state. The states with nonzero angular momenta have much higher energy jT^ ] and thus 
are irrelevant, therefore we will omit H below. Let us note now that the rest of the Hamiltonian (14) can be rewritten 
in the form 



E 



92 



dzjdz 



■\A{z,}\ + V{z,,z,}, 



where the potential energy V is given by 

v{z,,z.}^Y. E 



{Zi - Zk){Zi ~ Zj) 



4E 



92 



dzidz; 



log \zi 



El 



(15) 



(16) 



After redefinition of the wavefunction ^' |A{zi}|^' and rewriting Zi and Zi via real 2-vectors r^, we are back to the 
Hamiltonian defined in Eq. (Q) (apart from the constant term £'") with a — 1/2, m — 1/2 and n = I/tt. Evidently, 
the values of m and n give just the units of time and length in the problem, and can be rescaled arbitrary. The 
equilibrium probability distribution for the eigenvalues (for the density of eigenvalues equal to n) is: 



P{zi} ^ |*o{2;i}P const 



■ n 

j>k 



(17) 



and coincides (up to a constant normalizing factor) with the corresponding distribution for the general complex 
ensemble p8| , p^ (cf. also recent preprints ^9 20 1). Unfortunately it does not seem possible to establish an exact relation 
between the dynamics of general complex matrices and some Hamiltonian problem defined in terms of eigenvalues 
only; the problem is that a general complex matrix cannot be diagonalized by a unitary transformation and thus 
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angular variables cannot be eliminated from the dynamical problem |2l|. It is only the stationary distribution of 
its eigenvalues Pcompicx{zi} which coincides (up to a constant factor) with such a distribution for a normal matrix 
ensemble given by Eq. (^7|). However, the relations obtained above indicate that the correspondence between the 
KKLZ Hamiltonian and normal matrices is much deeper than just on the level of the stationary distribution. The 
supersymmetric structure of this Hamiltonian and the obtained relation with a relatively simple matrix model lead us 
to believe that the model under discussion might be a rather direct generalization of the Calogero-Sutherland model 
for a 2D case, and thus might even be amenable for the exact solution. However presently we are not aware of any 
exact way to treat normal matrix dynamics, and therefore will turn in the next Sections to some more elementary 
methods, which still make it possible to find some exact relations between the equal-time correlation functions of the 
KKLZ model at arbitrary value of the parameter a. 



IV. THE RELATION BETWEEN THE CURRENT-CURRENT AND DENSITY-DENSITY 

CORRELATION FUNCTIONS 



The SUSY structure of the Hamiltonian (Q) makes it possible to find, using the method used previously by Taniguchi, 
Shastry and Altshuler for the Calogero-Sutherland model (see also recent preprint |^^), an exact relation between 
the regular parts of the equal-time current-current lV°^(r) and density-density K{r) correlation functions defined as 

- r') = (J.(r) J^r'))"-^^; i^(r - r') = (p(r)p(r'))'-'=^. (18) 



In the definitions (|18| ) the singular (i.e. proportional to (5-functions) terms in the current and density correlators are 
excluded (i.e. r r'). It means, in particular, that the correlator K{r) is related to the inverse Fourrier-transform 
s(r) of the static structure factor 5(p) = / {dco) S {uj , p) as follows: 



K{r-r') ^ (0|^^(r-rj)^(r'-rfc)|0) = .s(r - r') - n5(r - r') 



(19) 



where n is the mean density of particles; the last term in ( [l9| ) is due to the fact that K{r) (contrary to s(r)) is the 
reducible correlation function). Similarly, the relation between the full current-current correlation function and its 
regular part is: 



Am? 



IE 



d 



5{v-vj) + 5{v~vj) 



d 



k,0 



dr, 



k,l3 



|0) 



Tlapiv - r') £;km<5Q/3(5(r - r') 



(20) 



where i^kin = ^ mean ground-state kinetic energy per particle (in D dimensions the last term in 

Eq. (|^) should be multiplied by 2/D). 

The trick to prove the relation between H™|(r) and K{v) is to substitute, following the standard definition of 
the current operator by the equivalent expression 



Ja(r) 



Then, using the commutation relations 



2m ^—^ 



4a 



(1 - 5,j)Dc,p{v, - Tj) - 5y ^ (I?a/3(r,; 



with Dap{r) ~ Vo-V^lnr, and identities (ji,a|0) = 0, we get 



K{r). 



(21) 



(22) 



(23) 



The relation ( |23| ) is applicable for any value of the coupling constant a (remember that a ^ 1 corresponds to the 
superfluid groimd state, whereas at very large a > acr ~ 35 the Wigner-type crystal is formed p^). One may 
wonder whether it is allowed to use the relation ( p3| ) at a > acr since in course of its derivation we have used space 
homogeneity of the correlation functions, which usually is not present in the crystalline state; however we will see 
below that the "crystalline" state present in our problem actually does not possess a true positional long-range order 
and thus does not break space homogeneity. 
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V. TRANSVERSE CURRENTS AND "INSTANTANEOUS" MAGNETIC RESPONSE 



We are going to obtain the momentum-space representations for the longitudinal and transverse current-current 
correlators defined as nQ^(q) = qaqfjll\\ (q) -I- {Sap — qaqp)Il±{q}. Performing the Fourrier transformation of the 
expression (^) and using the relation (^), we get 



n.Mq) = ^ / V^T^r" -^(P) + ^n'^ + -E^,Jap. (24) 

m-^ J {2ttY (p + q)"' q^ m 



d^p {p + q)aip + q)0 £.f^^^ , 2na^2lajl0 , 

r,2 



Here K{p) = >S'(p) — n is the Fourrier-transform of K{r) — . In order to obtain the transverse polarization function 
n_L(q) one has to contract Eq. ( p4| ) with the transverse projector Sap — qaQp'- 



2'Ka f cPp 



n±(q) = — / 



[(p + q)q]' 
(p + q)^9^ 



i^(p) -£;ki„. (25) 
m 



The q limit of the frequency-dependent nj^(q, w) is zero for any lj, leading to the same result for the equal-time 
correlator: n_L(q ^ 0) = 0. This condition makes it possible to determine -Ekin^ 

i?Mn = - ™ / yl^if (P) = -— [K{r = 0) - n^] . (26) 

Now, K{r = 0) = since the ground state (||) forbids 2 particles to be at the same point: 5'o(-.-, r, r, ...) = 0. 
Consequently, the expectation value of the kinetic energy in the ground state is given by 

„ ■nan uon 

i^kin = = ^, 27 

TO 4 



where wq = ^J^ng^^Jm is the plasma frequency. The relation (|2^) was checked numerically (see Appendix for the 
details) and found to be obeyed with a good accuracy for all values of a, including those there Quantum Hcxatic 
phase is formed, where the Galilei and rotational symmetries are both broken. We will discuss the implications of 
this result below, in Sec. VII. 

Then, Eq. (^ may be rewritten as 



^ , , 2'Ka f cPp 
n± q = 5" / -TTT^ 



[(p + q)q]' 1' 

(p + q)V 2 



(5(p) - n) . (28) 



If the orientational symmetry is unbroken (which is the case for not very large a), the correlation function S'(p) 
depends on the absolute value of the vector p only. Then the angular integral in Eq. (^8|) can be easily calculated, 
the result is: 

n±(g) = j\s(j>) -n)[l- ^) pdp. (29) 

At small momenta p the structure factor S{p) goes to zero: S{p — > 0) = np^ /2mujQ, therefore the low-momenta limit 
of the transverse correlator reduces to 

^^{q ^ 0) = ■ q\ 30 

OTO^ 

In the high-g limit the integral in Eq. ( p9| ) reduces to the one of Eq. (|2 

nx (z^«3 31 

TO^ 

In the special case a = 1/2 the current-current correlation function can be calculated explicitly. Indeed, the real-space 
density-density correlation function is known ||l^,|l^ : 

^1/2(0 ^n^[l~ exp(-7rnr2)] (32) 



(this expression was derived by direct integration of the total distribution function (|l^) over positions of all but two 
particles). The corresponding structure factor is 
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Si/2{p) = n[l - exp(-p^/47rn)], 



(33) 



which leads, via integration in Eq. (E9[), to 



n 



,1/2 



2m? 



1 - 



1 — exp(— (7^/47rri) 
q"^ /Ann 



(34) 



Note that the current-current correlations in the Calogero-Sutherland model were shown to be related to the cor- 
relations of eigenfunctions of the corresponding random matrices pO| ; we expect that a similar relation might exist 
between the correlation function ( p^ ) and some eigenfunctions correlations in the random ensemble of normal matrices. 

We have employed numerical simulations (cf. Appendix) to obtain the structure factor S{p) for several different 
values of a, and then used the obtained results to compute, using Eq. (p9|), the transverse current-current correlation 
function II±(p). The results are shown in Fig. 1. The peaks at p ^ 27m^/^ correspond to a roton minimum in the 
excitation spectrum that becomes deeper as a increases. 




FIG. 1. The structure factor S{p) (left axis) and the transverse current-current correlator n±(p) (right axis) computed 
numerically by the mapping of the Schrodinger equation with the Hamiltonian (^) onto a Langevin dynamics of classical 
particles (see Appendix). The curves correspond to a = 0.5, 1, 2, 5, 15 and 30. The increase of a leads to the increase of the 
peak height at p 27™^/^. 



VI. LONGITUDINAL CURRENTS AND LOW-Q PLASMON DISPERSION. 

Now we turn to the calculation of the longitudinal correlator n||(q). Contracting Eq. ( [2^ ) with the longitudinal 
projector and using (|2^), we obtain 

nil (q) - — n2 - nx (q) - — + — / {S{p) - n) [l ~ i^) pdp. (35) 

Alternatively, the value n||(0) may be found with the help of the continuity equation that leads to the relation 
q'^Il\\{q,uj) — uP'S{q,ijj) for the frequency-dependent correlators. At g ^ the whole spectral weight is saturated by 
the plasmon excitations with uj — ujq, so we get 
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9 n||(g) = / — ^ (36) 



leading to n||(0) = nuJo/'2,m in agreement with ( |3f: 

Eq. ( |35[ ) can be used to find the low-g dispersion of the plasmon energy to^q): using the continuity equation 
q'^Il\\{q,uj) — uo'^S{q,uj) together with the /-sum rule J{duj)LuS{q,Lu) = nq^ /2m, we convert Eq. ( |35| ) into the form 



q 



^Cj{q) = — ^ / {S{p) - n)v' dp, (37) 
oq^^ mn q^ Jq 



where 



j uj^S{q,uj)duj 
J ujJ(q, LUjduj 

As long as the spectrum ^(9) is nondecayable [i.e. Ima;(g) = 0], i^iq) coincides with Lu{q). 

Finally, using the relation S(q) = nq'^ /2mLL>{q) valid at g — s- and neglecting small S{q) compared to n in the r.h.s. 
of Eq. (|3^) , we arrive at 

uj{q) - Wo = -a- — . (39) 
4m 

Eq. (|3^) is very similar to the one derived |2^ ] in the lowest-order perturbation theory (over a) for the purely Coulomb 
2D Bose-liquid (the result of |Q contains extra numerical factor « 0.8 in the r.h.s. of Eq. (|3^)). Let us emphasize, 
however, that in the present model we are able to derive the relation (^ ) for any value of a, the only limitation was the 
low-g limit. In the field-theoretic language, the low-momentum excitation spectrum of the KKLZ Hamiltonian is given 
by the first-loop term in the diagram expansion. We believe that this remarkable property is due to the supersymmetry 
of the model, leading somehow to the cancellation of all higher-order diagrams. It is not yet clear to us whether the 
one-loop calculation is exact not only for the low-q spectrum, but for the spectrum and the correlation functions in 
general. To study this issue one needs to develop an explicitly supersymmetric second-quantized representation and 
the corresponding diagram technique for the KKLZ model. 

As long as a < 1, the quadratic dispersion uj{q) given by Eq. (^9|) produces a relatively weak minimum between 
ujo and the free-particle limit Wfroo('?) = q'^/2m. Such a spectrum is consistent with a totally superfluid (ris = n) 
ground state. With the increase of a the minimum in the spectrum becomes deeper; for a qualitative description of 
the spectrum one can use a kind of an interpolation formula 

u^iq) ^ - uJoaq^/2m+ {q'^/2mf. (40) 

The minimum of the r.h.s. of Eq. ( ^ ) touches zero at a = 2. It can be considered as an indication of possible 
spectrum instability and phase transition(s) in the large-a region. However the nature of Eq. ( ^o|) as interpolation 
formula of unknown (at large a) accuracy makes it impossible to use it for any exact conclusion about the behavior of 
the spectrum. In the next Section we will try to consider this problem starting from the opposite limit of extremely 
large a. 



VII. "CRYSTALLINE" GROUND STATE: QUANTUM HEXATIC WITHOUT SHEAR MODULUS 

At very large values of a > Ofcr ~ 35 the ground state is expected to possess symmetry breaking realized in 
the formation of triangular lattice |p3|| . The static structure factor S{p) can be determined knowing the ground 
state wavefunction only. Its calculation involves integration over all particle's coordinates with the weight given 
by ^Polri}, which is formally equivalent to the calculation of the structure factor of a classical 2D crystal with a 
logarithmic repulsion between particles V{r) = 4q; • Inr and "temperature" T = 1 (cf. Appendix). Such a classical 
crystal has a finite shear modulus fXd ^ an and is expected to undergo dislocation-mediated melting when a decreases 
below some critical value ai. At a > ai it is stable with respect to dislocation formation, but does not possess long- 
range translational order (as any other 2D crystal at finite temperature). It means that the structure factor S{p) does 
not have Bragg peaks near the principal inverse lattice vectors Pa, but shows only power- low singularities: 

S{p)^\p-Par' (41) 
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with rj . Let us now go back to our original problem, which is a quantum 2D crystal. Suppose that it has finite 

shear modulus /i (now without subscript 'cl'), then the spectrum of the transverse sound would be linear, as usual: 
'^(9) ^ Vm*? where q is the wavevector of a shear deformation on a top of the crystal state. (In terms of momenta 
p of the original Bose particles it corresponds to ^^^(p) ~ |p — Pal for any a). Now we can calculate the zero-point 
fluctuations of lattice displacements {v?) ~ / d?quj^^{q) and find it finite, which means a finite Debye- Waller factor 
and thus the presence of Bragg peaks in the structure factor S(j)) - in contradiction with the above-obtained Eq. (|4l|). 
The only way out of this contradiction is to suppose that the quantum solid with our very special interaction given 
by Eq. (^) has actually zero shear modulus: ^ = 0. More precisely, we need to suppose that uj{q) cx q"^, which leads 
then to a logarithmically divergent (u^), in agreement with Eq. (pl|). It means that the "crystalline" ground state 
expected to exist at a > ai should be more adequately called "Quantum Hcxatic" (QHX), which is (in close analogy 
to the classical hexatic phase introduced by Halperin and Nelson within their picture of dislocation-mediated thermal 
melting of 2D crystals Q) the phase without shear modulus but still with anisotropy of the excitation spectrum: 

Wi(p)^(p-p„)VM (42) 

(M is some parameter with dimensionality of mass), whereas in the isotropic liquid excitation spectrum depends on 
the momentum absolute value |p| only. 

The absence of a shear modulus in the high-a ground-state of the Hamiltonian (^ can also be shown by direct 
calculation which can be conveniently performed representing potential energy of interaction between particles as 

^KKLz = — E(E(ir^j +— E^(r^-.) (43) 

Now one can calculate the variation of the potential energy (actually the second sum in (^ ) does not contribute) 
corresponding to an infinitesimal uniform shear deformation eap of the regular lattice. Then one finds that it does not 
contain the term proportional to e^^, i.e. the expansion starts from the term (Veq/j)^, which means that /i(<?) oc q^. 

The absence of a genuine crystalline state in the KKLZ model at arbitrary strong interaction a — *■ cxd stems from 
its very special form and certainly is not present in a generic 2D Bose system. However, in a more realistic system 
QHX may exist as an intermediate ground-state realized at moderately large repulsion. 

Let us consider now the response of the QHX state to a magnetic field h. Since particles in this state are almost 
localized, we can neglect (in a first approximation) the multi-particle exchange contributions and estimate suscep- 
tibility as X ~ —{v?) ~ — \n[nL?) (where L is the system's size). Comparing this expression with the magnetic 
susceptibility of a superfluid, which is proportional to UgL^, one concludes that = in the QHX state. 

Now we need to address some apparent paradox: i) according to Eq. (|lO|), normal density n„ = n — Us is proportional 
to the integral n^(g — > 0,uj)dLj/iLj where the expression under the integral is positive; ii) the equal-time current- 
current correlation function Il±{q) = /q°° Ilj_{q,uj)du!/Tr vanishes at g — ^ (cf. Eq. (|2^)); iii) then, how it is possible to 
get nonzero n„ (and even cancel Ug completely)? In fact, to get nonzero n„ we need to allow for the singular behavior 
of the function Il±(q,uj) at low q,uj (in other words, the limiting procedure q —> and integration in Eq. ( p^ should 
not commute with each other). As a simple example one can consider Il{q,uj) oc f{q^ jmuj) with f{x) = x/{\ + x^). 
However, such a behavior of H((7, lo) is only possible if there are some gapless excitations in the system, which is just 
the case in the QHX state (cf. Eq. (^)). In general, an appearance of gapless excitations with nonzero momentum 
(low-momenta excitations in our system always have large plasmon gap loq) may only be a result of a phase transition, 
which is associated with spontaneous breakdown of the Galilei invariance. Generally speaking, the Galilei invariance 
may be spontaneously broken at a = ac < which would mean the existence of two phase transitions and some 
intermediate state between the SF and the QHX states. 



VIII. BETWEEN SUPERFLUIDITY AND QUANTUM HEXATIC: SPECULATIONS 

We found in previous Sections that at small a the ground-state possesses total superfluidity {ug = n) and the 
excitation spectrum beyond the ground state is gapful (cf. Eq. (^o|)). At very large a > acr ~ 35 the ground state 
is not superfluid at all, Ug — 0, and the excitation spectrum is very soft (cf. Eq. (^^). How does the transition(s) 
between these two qualitatively different ground states occur? The most trivial possibility would be that of the direct 
flrst-order phase transition at a = aa-. This would also mean that the flrst-order melting/crystallization transition 
happens in the classical 2D Coulomb gas at the interaction parameter F = 4a « 140. However, the persistence of 
total superfluidity up to a « 35 looks rather improbable since at so large values of a the characteristic energy of the 
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interaction between particles, (7^/27r, exceeds the kinetic energy i5kin = Tran/m by more than 2 orders of magnitude. 
Note also, that the above-cited value of E'kin is about 100 times larger than the naive estimate n/m; it means that 
particles are almost localized in rather small regions of space (much smaller than the mean interparticle distance 
n~^/^) and diffuse from these regions on much longer time scale only. Therefore we suspect that at least two phase 
transitions take place between the totally SF state and the QHX state. 

As we discussed above, the QHX quantum ground-state is seen as a 2D Berezinsky crystal with logarithmic inter- 
action as long as only the equal-time density correlations are considered. The nature of melting in such a 2D classical 
crystal is not quite clear. The old numerical simulations reported weak first-order transition p^, but the recent 
simulations on the 2D systems with other types of interactions reveal that to discern between 1st and 2nd order 
melting one probably would need to simulate systems with a huge number of particles, which was not yet done for 
the system of interest here. 

If melting of such a 2D crystal occurs via Halperin-Nelson mechanism, then there exists a usual hexatic state which 
possesses quasi-long-range orientational order, i.e. power-law decay of the orientational order parameter 

(exp[6«(6'(0) - e{Y))]) (X r-'T' (44) 

without any translational order. The corresponding quantum ground-state would be characterized, apart from orien- 
tational quasi-LRO, by the isotropic structure factor S{p) and by the existence of some gapless mode related with a 
slow rotation of the 9{y) angle. The other possibility is that the melting of the 2D crystal with log-interaction occurs 
via a first-order phase transition, without hexatic phase; even in this case there is a room for a phase transition related 
to the breakdown of the Galilei invariance, which would be manifested by an appearance of the gapless mode in the 
excitation spectrum. Indeed, the absence of normal density n„ in the ground-state of the system with purely potential 
interactions is the consequence of the Galilei invariance of the Hamiltonian, which is certainly respected in the low-a 
limit. An appearance of non-zero n„ is thus related to spontaneous breakdown of the Galilei invariance, which may 
occur at some value of a which we denote as ac- As it was explained above, we expect ac < ctcr, although we still 
cannot prove this inequality. 

To conclude this Section: if there exists any additional ground-state sandwiched between SF and QHX, then the 
spectrum of excitations beyond such a state contains a gapless mode with isotropic spectrum: uj{p) cx [p — pa)^. Such 
a mode may be orthogonal to the density mode and, thus, it may be invisible in the usual structure factor S{q^uj). 



IX. CONCLUSIONS 



We have studied the model Hamiltonian of the two-dimensional strongly interacting Bose-liquid, which can be 
considered as a natural generalization of the Calogero-Sutherland Hamiltonian for the 2D case. This Hamiltonian is 
shown to possess the non-relativistic supersymmetry. It makes possible to derive a number of exact relations between 
the equal-time density-density and current-current correlation functions, and to find exactly the low- momentum 
excitation spectrum of the model, which is given by the one- loop correction to the mean-field result. The ground 
state is shown to exist in at least two qualitatively different phases as function of the interaction parameter a: the 
superfluid state at small a and the quantum hexatic state at very large a. The arguments are presented in favor of 
the existence of a third state, intermediate between the two above-mentioned ones. For the special case of a — 1/2 the 
considered model is shown to be equivalent to the quantum mechanics of normal complex matrices with a Gaussian 
weight, exactly in the same way as the Calogero-Sutherland model at the special values of its coupling constant is 
equivalent to the dynamics of the Wigner-Dyson matrix ensembles. 

The model we have studied seems to be very interesting as a natural approximation for the problem of 2D Bose- 
liquid with logarithmic pair-wise interaction, which is widely considered in relation with the composite-boson approach 
to the description of the Fractional Quantum Hall state |8[|2^, as well as in relation with statistics of vortex lines 
in HTSC superconductors The key physical problem here is the problem of existence of non-superfluid liquid 

ground-state. We have shown within the considered model the existence of at least one nontrivial ground-state, which 
is quantum hexatic, i.e. anisotropic non-superfluid liquid. It is still an open question if the destruction of superfluidity 
may occur without any explicit symmetry-breaking in the real space, i.e. if it might be that the only spontaneously 
broken symmetry is the Galilei invariance. To answer this question it is necessary to study the spectrum of lowest- 
energy excitations at large wavevectors q ~ l-K^/n and at large a € [2, 35] and look for the appearance of some gapless 
mode; note that such a mode, if it exists, may be orthogonal to the density-oscillation mode, and thus may be absent 
from the dynamic structure factor S{q,Lu). 

The exact results obtained make it possible to suspect that the considered model can actually be exactly solvable, 
at least at some special values of a. There are several directions in which one could try to extend and develop the 
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results presented here: i) to find the time-dependent correlation functions, at least for a = 1/2, by the solution of 
the normal-matrix problem; ii) to generalize the correspondence between our Hamiltonian and the normal matrix 
problem for other values of a; iii) to develop an explicitly supersymmetric second-quantized many-body theory and 
the corresponding diagram technique for our Hamiltonian; this might be very helpful in order to generalize our result 
for the spectrum, Eq. (|39|), for an arbitrary momenta q. A completely different direction of study of this model is 
provided by its equivalence to a Fokker-Plank equation of multi-particle stochastic dynamics (cf. Appendix), which 
makes it possible to use c/assica/ Monte-Carlo simulations in order to study some properties of the quantum many-body 
problem in question. 

We are grateful to V.B.Geshkenbein, A.Gorsky, A.I.Larkin, A.Yu.Morozov, A.M.Polyakov and A.S.Schwartz for 
useful discussions. We are especially indebted to L.B.loffe for his collaboration on earlier stages of this research and 
many useful advises. The hospitality of the Isaac Newton Institute where this study was initiated during the visit 
of one of us (M.V.F., August 1994) is gratefully acknowledged. This research was supported by the INTAS-RFBR 
grant # 95-0302, Swiss National Science Foundation collaboration grant # 7SUP J048531, and DGA grant # 94-1189 
(M.V.F.). 



APPENDIX: SUSY, LANGEVIN DYNAMICS AND NUMERICAL SIMULATIONS 

The relation between supersymmetric quantum mechanics and Fokker-Planck equation of stochastic classical dy- 
namics is well-known now |l^,^ (cf. also recent preprint ||2^). In this Appendix we will show how to use that 
equivalence in order to implement fast numerical simulation of the many-body quantum-mechanical problem we are 
interested in. 

Consider a system of iV coupled Langevin equations: 
with a Gaussian distribution of random forces C(t): 



iiA^)iK(m = 2-iT5,k5c.p5{t - t')- (46) 

Here T is temperature of a thermal bath, 7 is a damping coefficient which can be set to 1 by a proper rescaling of 
time. 

The Fokker-Planck equation for the probability density function, 

(47) 

after the substitution V{ri} — e^^/^^\I'{ri} reduces to an imaginary-time Shrodinger equation 

1 d^W 1 / dW 



0^ \ 

h\ 



* ) . (48) 



The many-particle ground state of Eq. ([48| ) has zero energy, which is a consequence of a supersymmetry |13| ], and 
wavefunction = const e"'^'^'"^/^^. 

An equivalence of Langevin, Fokker-Planck and Shrodinger approaches implies that the averaging over many-particle 
ground state ^"0 is equivalent to the averaging over the Langevin dynamics of particles (with respect to random force 
and initial particles distribution) p^: 



vI/o|F{rJ|*o) - (i^{rj)^. (49) 



The relation (49) is very important as it makes possible to compute ground-state averages by numerical simulation 
of Langevin dynamics with a proper interparticle interaction W which is technically substantially simpler than direct 
quantum Monte-Carlo simulations. The problem at hand, given by the wavefunction (^, corresponds to a stochastical 
classical dynamics of a one-component 2D Coulomb gas governed by the following potential energy: 

W{vi\ = -4aT^ln \yj -Vk\+ iT^oTn'^rl. (50) 

j>k i 
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We have performed Monte-Carlo computer simulations of the Brownian dynamics (^) with the potential for 
the systems oi N = 109 particles with different values of the parameter a. The results obtained for the structure 
factor S{p) are shown in Fig. 1. Also shown are the curves for the correlation function of transverse currents, Il±{p), 
calculated from S{p) according to Eq. (|29|). 

We have also used this simulation in order to check Eq. ( p7| ) for different values a. Note that the general relation 
( p9| ) can not be applied to compute the mean value of the kinetic energy since the operator of the kinetic energy 
contains space derivatives that have no immediate meaning in the classical Langevin dynamics. Fortunately, instead 
of i?kin one can compute the mean ground-state potential energy E'pot using the fact that the ground-state energy is 
zero. Epot may be easily obtained by Langevin simulations using Eq. since the operator of the potential energy 

t^KKLZ + + ^"^^^^ ''f contains particles coordinates only. 

Numerical simulation have showed that Eq. holds for all values of a from < 1, that corresponds to the SF 
state, to 100, corresponding to the CR state, within the precision of computation. 
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